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ABSTRACT 

Dust extinction is generally the least tractable systematic uncertainty in astron- 
omy, and particularly in supernova science. Often in the past, studies have used the 
equivalent width of Na I D absorption measured from low-resolution spectra as proxies 
for extinction, based on tentative correlations that were drawn from limited data sets. 
We show here, based on 443 low-resolution spectra of 172 Type la supernovae for 
which we have measured the dust extinction as well as the equivalent width of Nal D, 
that the two barely correlate. We briefly examine the causes for this large scatter that 
effectively prevents one from inferring extinction using this method. 



1 INTRODUCTION 

The ubiquity of dust in galactic environments, most notably 
in star forming regions where many supernovae (SNe) ex- 
plode, makes extinction correction one of the most pervasive 
and probably the least tractable systematic uncertainty in 
the study of SNe. When an object's intrinsic colors are un- 
known, as occurs frequently with rare and new classes, one 
is bound to look for indirect measurements of the amount of 
extinction and reddening that need to be taken into account. 

The absorption doublet of sodium NalD, at 5890 and 
5896 A, is a well-known tracer of gas, and its strength is 
generally expected to indicate the amount of dust along the 



cesses measured for these stars, with a noticeable scatter. 



line of sight. Ferlet et al. ( 1985 1, for example, show that the 



column densities of Na I D correlate with those of hydrogen 
in the diffuse interstellar medium. However, line saturation 
often makes density estimations arduous (e.g., ,Mugglestone| 
|fc 0'Mara|1965[ ). Variations in dust-to-gas ratios in different 



galaxy types (see, for example, Issa et aL]|1990 Lisenfeld & 



Ferrara 1998 1 make the jump from sodium to dust some- 



what more uncertain. Depletion of metals on dust grains 
(e.g.. Savage & Mathis 19791 further complicates the pic- 



ture. It would therefore be quite surprising to find a tight 
correlation between the strength of sodium absorption and 
dust extinction. 

Richmond et al. ( 1994 1 have shown, using 57 high- 



resolution stellar spectra compiled from the literature, that 
the equivalent width (hereafter EW) of the individual com- 
ponents of NalD does indeed correlate with the color ex- 



Munari & Zwitter ( 1997 1 add a body of 32 stars observed 
with a somewhat lower resolution {R ~ 16,500 vs. R ~ 
60,000-600,000); they constrain the nonlinearity of the rela- 
tion, extensively discuss its pitfalls (such as when there are 
multiple absorption components, or when the extinction is 
higher than E{B — V) ~ 0.4 mag), and find that the preci- 
sion is limited to about SE{B — 1/) « 0.05-0.15 mag. They 
also show how at EW > 0.5 A the NalDi line saturates and 
the relation fiattens. 

Type la SNe — understood to be the result of the ther- 
monuclear disruption of a white dwarf (see, for example, 
Hillebrandt fc Niemeyer|2000 for a review) — are very effec- 
tive distance indicators through "standardization" of their 
light curves; they were used to discover the accelerating ex- 
pansion of the Universe (Ri ess et al.|1998 Perlmutter et al. 
l999|and to accurately measure the Hubble constant ( [Riess 
et al.|20lT and references therein). Since their luminosity is 



found to correlate well with their light-curve shape ( Phillips 
|1993[ ), various methods fit those light curves to templates in 
order to extract a distance (e.g., Riess et al.|1999 Guy et al 
|2005' 



Conley et al. 20081. In the process, dust extinction 
effectively a nuisance parameter — is treated in various 
ways, either fit separately or as part of a more generic color 
dispersion term. 

Barbon et al.l (19901 and iTuratto et al.l (|2003|) used 



this property of SNe la to compare the EW of NalD mea- 
sured in low- to medium-resolution spectra to the extinction 
they derived from light-curve fitting. This would extend the 
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method to the significantly more common occurrence, when 
one does not have high-resolution spectra and the doublet is 
blended. Using few spectra in each case (6 and ~30, respec- 
tively), they derived tentative scaling relations which have 
been widely used in the literature over the last two decades. 
However, as we show below through a similar analysis of a 
much larger sample, these relations are of little predictive 
power. While a correlation indeed exists between extinction 
and the EW of NalD, the very large scatter makes it es- 
sentially useless when only low-resolution measurements are 
availablaM 



2 DATA AND MEASUREMENT 

For over a decade the Lick Observatory SN search (LOSS; 
Li et al.|2000l|Filippenko et al.|26oT| with the Katzman Au- 
tomatic Imaging Telescope (KAIT) has been one the most 
successful systematic discovery engines of nearby SNe. KAIT 
has also conducted follow-up photometry of hundreds of 
SNe, along with the Nickel 1-m telescope at the Lick Ob- 
servatory. Moreover, the LOSS team has obtained spectra 
of these SNe with multiple telescopes, but mostly with the 
Kast double spectrograph (Miller & Stone 1993 1 mounted on 



has been widely tested and vetted, far beyond our needs for 
this work. 



2.2 Na I D Equivalent Width 

From our database of 1613 spectra of 743 SNe la (Silverman 
et al., in preparation), we choose every spectrum of SNe 
having an MLCS estimate, resulting in 724 spectra of 239 
individual SNe. In order to simplify the analysis we further 
cull the sample with the following criteria. (1) Objects that 
are found to be normal (i.e., non-peculiar) SNe la based on 



the Lick Observatory 3 m Shane reflector, the Low Resolu- 
tion Imaging Spectrometer (LRIS; Oke et al.|1995 I mounted 
on the 10 m Keck telescopes, and the Deep Imaging Multi- 



Object Spectrograph (DEIMOS; Faber et al. 20031 on the 



10 m Keck II telescope. The photometric data are presented 
by Ganeshalingam et al. (120101), and the spectroscopy by 



Silverman et al. (in preparation). In addition, we have com- 
piled photometry from the Calan/Tololo SN search (Hamuy 
et al.|1996 \ and the Harvard Cf A samples ( Riess et al.|1999 
Jha et al.|2006l|Hicken et al.||2009[ ). 



2.1 Extinction Estimation 



matching with the SN identification code (SNID; Blondin 
& Tonry 20071. (2) Only spectra obtained with the three 
main spectrographs named above. (3) Redshift z > 0.005, 
in order to avoid contamination of the EW measurement 
with observer-frame absorption from our own Galaxy. (4) 
MLCS goodness of fit of < 2, and estimated uncertainty 
in Av of 5Av < 0.2 mag. These cuts result in a sample of 443 
spectra of 172 SNe. This sample is still more than an order 
of magnitude greater than in previous studies. Our results 
below are independent of the precise cuts applied (such as 
the redshift limit or minimal goodness of fit). 

The procedure for deriving the EW of Na I D for every 
spectrum is as follows. The spectrum is inspected visually, 
and if a line is noticeable its edges are manually marked. We 
perform numerical integration of the line over this marked 
range, dividing the result by the integration over the pseudo- 
continuum which is assumed to be linear. We divide the 
spectrum by a low-order polynomial fit to the rest-frame 
range of 5800-6000 A and measure the noise of the spectrum, 
A*', defined as the 5a clipped standard deviation around the 
flattened spectrum. The uncertainty in the EW is assumed 
to be dominated by A'^. We calculate the uncertainty, SEW, 
as SEW = Nx W/2, where W is the width (in A) of the line 
(the factor of 2 accounts for the approximately triangular 
shape of the feature). In cases where no line was detectable, 



Ganeshalingam et al. (in preparation; see also |Gane-| r ,i_ i-' ■ , /or, X 
: 2 11 . \ f t- 1 I _ — I width for the hue m our spectra (20 A 



shalingam et al. 2010 1 fit every SN la to templates using 



ty \ 

h 



the multicolor light-curve shape fitter, MLCS2k2.v006 (here- 
after MLCS; [Riess et al.|199"6{|Jha et al.|2007| ), providing an 
estimate of the amount of extinction suffered by the SN in 
its host galaxy. Milky Way reddening is removed by using 



the dust maps of Schlegel et al. ( 1998 1. Two independent fits 



are performed, using a galactic line of sight prior extinction 
law with the Milky Way average of Rv = 3.1 or a steeper 
Rv = 1.7 (e.g., |Wang et al.][2ro9l [Hicken et aL][2009l and 



references therein). While the two sets give systematically 
different values for the amount of extinction Av, the results 
are highly correlated and do not affect our conclusions be- 
low. For simplicity, we discuss only the values derived with 
Rv = 3.1, which is consistent the value of Rv = 2.8 ± 0.3 



found by Chotard et al. (2011 1 



While there may be some debate as to the accuracy of 
the derived extinction values with MLCS due to the depen- 
dence on the prior dust law and to the degeneracy between 
intrinsic color variance and external reddening, the precision 



2.3 Variability 

For 85 of the SNe in our sample we have more than one 
spectrum, the majority with three or more spectra. For ev- 
ery such SN, we examine the stability of the EW of Nal D by 
comparing the result for a given spectrum to the weighted 
mean for all spectra of that object. In only four cases (out of 
254 spectra) do we find an EW that is inconsistent by more 
than 2a with other spectra of the same SN. In two of those 
cases we find that this is due to a slight underestimation 
of the uncertainty. For the remaining two SNe, SN 2006bq 
and SN 2006cm, the variation in EW seems to arise from a 
varying amount of host-galaxy light spilling into the spectro- 
scopic slit. This finding raises some concerns regarding the 
reliability of the EW of Na I D as a good proxy for dust ex- 
tinction, as measurements for the same SN can, on occasion, 
vary significantly. 



^ See also a similar discussion by N. Elias de la Rosa at 
http: //online .kitp.ucsb . edu/online/snovae07/eliasdelarosa/ 



and 



Elias- Rosa 



120071. 



^ A table including all of our measurements is available at 
http: // astro .berkeley. edu/~dovi/Poznaiiski2011_SNjnad. dat 
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Figure 1. In red, extinction along the line of sight, Ay, from template fitting to the light curves vs. the equivalent width of NalD as 
measured from low-resolution spectra of SNe la. The large blue squares are weighted averages in bins of equal sample size. Nondetections 
are marked by arrows located at the la limit. In addition we show the various correlations proposed by |Barbon et alT] ( |1990[ l, |Richmond| 
|et al.| ( |1994^ , and [Turatto et al.| | |2003[ two different lines) . While there is indeed a correlation as seen in the mostly monotonic rise of 
the binned data, the scatter is far too great to make this spectral feature an effective predictor of dust extinction. 



3 CORRELATION? 



4 CONCLUSIONS 



In Figure [T] we show in red the value of Ay derived from 
MLCS compared to the EW of NalD. Nondetections are 
marked with arrows at the la limit. We bin the EW mea- 
surements, with a variable bin size in order to have the same 
number of spectra per bin (about 50). We plot with large 
blue squares the weighted mean (in EW and Ay) for each 
bin. Here a shadow of the correlation arises, but the disper- 
sion around it is too great to make it a useful tool. Our mea- 
surements are effectively consistent with all the tentative 
correlations sugg ested by [Barbon et ah] ( |1990| ), |Richmoiid| 
et al. (19941, and Turatto et al. (20031, as shown in the fig- 
ure. Our best fit correlation is S(_B-y) = 0.43 EW[A]-0.08, 
with a systematic scatter of 0.3 mag (1(t). A similar scat- 
ter is observed when plotting EW vs. MLCS Ay values for 



7?v = 1.7 or the SALT2 c term (|Guy et al.||2007l 



Our finding is that the EW of NaT D, as measured from low- 
resolution spectra, only very weakly tracks the amount of ex- 
tinction suffered by a SN, at least as long as the extinction is 
in the range we probed. Ay < 3 mag. Many effects conspire 
to minimize the amount of inference one could make. SN 
spectra are typically contaminated with host-galaxy light, 
and the continuum emission will dilute the line strength. 
This contamination is highly dependent on the slit width 
and position, as well as on observing conditions. 

In addition, low-resolution spectra cannot resolve the 
NalD doublet and therefore the different curves of growth 
of the two components are blended. Variations in gas-to- 
dust ratios in particular regions or galaxy types, multiple 
gas clouds at different velocities along the line of sight, and 
line saturation at large optical depths, all contribute to the 
theoretical and practical explanation for our results. 



In Figures [2] and [3] we illustrate the two most extreme 
outliers in our sample, SNe 1995E and 2006et. Both are well 
fit by MLCS, with Ay of 2.4 and 0.5 mag (respectively), but 
their Na I D EW values are about 0.7 A and 6 A. While these 
are far from typical. Figure [T] shows that they are clearly 
not unique, and they demonstrate the risk associated with 
deriving extinctions from NalD EW measurements in low- 
resolution spectra. 
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Figure 2. Left: Multicolor light curves of SN 1995E together with the best-fitting MLCS templates, from which a large extinction of 
Ay = 2.4 mag is derived. Right: Low-resolution spectrum of SN 1995E (inset shows the NalD feature in detail; blue lines mark the 
wavelengths of the doublet), indicating that the EW of NalD is rather small, about 0.7 A. In red we show what the line should look like 
in order to have the EW predicted by the relation of [Richmond et al.| | |1994] |. 




Figure 3. Left: Multicolor light curve of SN2006et together with the best-fitting MLCS templates, from which a moderate extinction 
of Ay = 0.5 mag is derived. Right: Low-resolution spectrum of SN2006et (inset shows the NalD feature in detail; blue lines mark the 
wavelengths of the doublet), indicating that the EW of NalD is extremely large, about 6 A. In red we show what the line should look 
like in order to have the EW predicted by the relation of [Richmond et al.| | |199^ . 



